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ARTICLE INFO ABSTRACT 








Keywords: Treatment in early intervention services (EIS) seems superior to treatment as usual on several outcomes, but the 
Psychosis extent of heterogeneity in response is unclear. In this study, treatment response trajectories up to 2 years in first- 
Schizophrenia episode psychosis (FEP) patients enrolled in an Italian early intervention service (EIS) have been quantified. The 
ies model 24-item Brief Psychiatric Rating Scale (BPRS) was used to quantify treatment response up to 2 years in 129 


participants. Conditional growth modeling and latent class growth analysis were used to test changes over time 
in the BPRS and separation into independent classes over time. Group differences were tested on socio-demo- 
graphic and clinical variables known to be related to outcome in psychosis. Scores on the BPRS showed a 
statistically significant decrease in overall scores across all tested models. Four trajectories were identified across 
2 years. Most patients showed a progressive decrease in the BPRS scores; a scant fraction showed a more stepped 
decrease from very high levels of psychopathology. No potential predictor was statistically related to the time 
course of BPRS scores. Most patients that undergo treatment within an EIS are characterized by amelioration, but 


Early intervention 


patients that have higher baseline scores of psychopathology require more intensive treatment. 





1. Introduction 


Schizophrenia spectrum disorders import a high burden of care, 
reduced quality of life, enhanced death rates, and elevated personal, 
family and societal costs (Awad and Voruganti, 2008; Laursen et al., 
2014; Cloutier et al., 2016; Jin and Mosweu, 2017; Watson et al., 2018). 
There is evidence that a combination of drug and psychosocial treat- 
ment may improve the long-term outcome of patients with schizo- 
phrenia spectrum disorders (Harding et al., 1987; Volavka and 
Vevera, 2018). There is also some evidence that early course patterns 
within the first 2 years of treatment are the strongest predictor of 
medium- (15 years) and long-term (25 years) outcome (Harrison et al., 
2001; Lally et al., 2017). In particular, early treatment response was 
related to long-term outcome (Emsley et al., 2008), and, in general, the 
duration of untreated psychosis (DUP), i.e. the time from manifestation 
of the first psychotic symptom to initiation of adequate treatment, was a 
predictor of outcome, with longer DUP related to poor general symp- 
tomatic outcome, poor social functioning and lesser likelihood of re- 
mission (Penttilä et al., 2014). However, the links between DUP and 
long-term outcome might be a reflection of the better prognosis of cases 
with acute onset when compared to cases with insidious onset 


(Volavka and Vevera, 2018). In first-episode psychosis (FEP), predictors 
of relapse following discontinuation of treatment after remission from 
the episode were male sex, prior psychiatric admission, and premorbid 
adjustment and functioning (Bowtell et al., 2018). Age of onset has a 
small but statistically significant impact on the outcome of patients 
diagnosed with schizophrenia, with younger age of onset being related 
to more hospitalization and poorer global outcome (Immonen et al., 
2017). However, age of onset might be confounded by gender differ- 
ences, with men having younger age of onset than women (Cascio et al., 
2012), and more likely than women to suffer from comorbidities that 
impact on outcome (Riecher-Réssler et al., 2018). Men, for example, are 
more likely than women to develop alcohol and substance use disorders 
(SUDs), which occur between 20% and 65% of samples depending on 
the setting and the type of SUD, with SUDs consistently related to 
poorer symptomatic outcome and greater chance of hospital admission 
as well as greater risk of violence, victimization and suicide 
(Kerfoot et al., 2011; Bennett et al., 2017). Family history of psychiatric 
disorders is a known risk factor for schizophrenia, however it is weakly 
related to global outcome (Käkelä et al., 2017). Nevertheless, when the 
investigation is restricted to a family history of psychosis, a relatively 
small but statistically significant association with the long-term 
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occupational and global outcome is found in patients with schizo- 
phrenia (Käkelä et al., 2014). Besides age, sex, family history of psy- 
chosis, DUP, premorbid functioning, and SUDs, dropout of treatment 
was related to outcome in patients diagnosed with schizophrenia- 
spectrum disorders. Dropout was consistently associated with a reduced 
chance of improvement, and the severity of the symptoms before the 
last visit was related to a greater chance of dropout (Rabinowitz and 
Davidov, 2008). 

Since the 19th century, clinicians agreed that treatment would be 
the most effective at the early stage of the mental illness. Among others, 
influential authors, such as Winslow (1860), who was specialized in 
forensic psychiatry (Winslow, 1858), and Jean-Etienne Dominique 
Esquirol (1845), one of the fathers of modern psychiatry 
(Charland, 2015), held that mental disorders are preceded by a period 
of incubation during which early symptoms, including changes of ha- 
bits, are prone to be overlooked. These authors suggested treating these 
early symptoms as soon as possible. However, during the 19th century 
patients were rarely treated in a timely fashion, also because of the lack 
of effective treatment, as summarized in a recent review on the topic 
(Chau et al., 2018). In modern times, converging evidence suggests that 
early intervention in psychosis may have a favorable impact on re- 
mission of positive and negative symptoms and functional recovery 
(McGorry, 2015; Santesteban-Echarri et al., 2017). Across 10 rando- 
mized clinical trials involving 2176 patients with FEP, treatment in 
early intervention services was superior to treatment as usual in total 
symptom severity (standardized mean difference [SMD], —0.32; 95% 
CI: —0.47 to — 0.17), positive symptom severity (SMD: — 0.22; —0.32 
to —0.11), negative symptom severity (SMD: — 0.28; —0.42 to —0.14), 
and in school or work functioning (risk ratio: 1.13; 95% CI: 1.03-1.24) 
(Correll et al., 2018). 

Most studies insofar had investigated the outcome of the first epi- 
sode of psychosis assuming a linear trajectory from baseline to the 
considered endpoints. However, there is increasing evidence that im- 
provement in patients with FEP is heterogeneous, with patients with 
better premorbid functioning and higher cognitive scores showing the 
most improvement (Levine et al., 2010; Levine and Rabinowitz, 2010), 
and those with most severe negative symptoms at inception having a 
reduced chance of improvement (Horton et al., 2015). Overall, those 
who continue to improve are more likely to stay in treatment compared 
to those who only experience temporary relief, who are more likely to 
drop out of treatment (Kong and Chen, 2016). Differences in the tra- 
jectories of improvement in FEP patients may be related to predictors 
that are amenable to intervention, thus their investigation may offer 
clues to improve the course of these disorders. 


1.1. Aims 


This study was set out to investigate the outcome of FEP patients 
who were treated for 2 years in an early intervention service. We aimed 
at a) establishing the impact of known predictors of outcome in schi- 
zophrenia spectrum disorders; and, b) estimating empirically-derived 
trajectories in change over time of the severity of symptoms to quantify 
the extent of heterogeneity in treatment response and to identify sub- 
group characteristics in treatment response. Mixed models of growth 
curves and latent class growth analysis were applied to pursue these 
aims of the study. 


2. Methods 


The study was carried out in the Programma2000 service, the early 
intervention service operating under the Health Authority of the ASST 
Niguarda Hospital of Milan (Cocchi et al., 2008). Data were collected 
during the routine assessment of the patients participating in the Pro- 
gramma2000 and, more specifically, of those enrolled from June 1999 
to December 2014. The study protocol complies with the guidelines of 
the 1995 Declaration of Helsinki, as revised in Tokyo in 2004, and 
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further revised in Fortaleza, Brazil, in 2013 (World Medical 
Association, 2013). All included patients gave their informed consent. 


2.1. Participants 


Programma2000 is a multi-modal, community-based program on 
early intervention in psychosis settled in Italy, Milan, since 1999. 
Programma2000 operates as an assessment and treatment hub for 
several referral agencies, ranging from institutionally mediated path- 
ways (e.g. primary care, district MH, school counseling, emergency 
rooms) to spontaneous help-seeking individuals (e.g. self-referrals or 
family referrals in response to awareness campaigns). The served 
catchment area includes approximately 200,000 inhabitants. Criteria 
for referral and preliminary evaluation are: age between 17 years and 
30 years and help-seeking for early or impending psychosis. 

Participants who were included in this study were patients diag- 
nosed with FEP, operatively defined as: a diagnosis of schizophrenia or 
related syndromes (such as F20-29 in ICD-10) according to the ICD-10 
(World Health Organization, 1992), with duration of untreated psy- 
chosis (DUP) less than 24 months, which is considered the limit to start 
an effective early intervention protocol of care (Birchwood, 2000). The 
majority of this group includes subjects referred to Programma2000 
immediately after a first contact with any public mental health service 
of the catchment area and contextually diagnosed as suffering from 
non-affective psychosis prior to admission. 

Exclusion criteria were: comorbid medical or neurological dis- 
orders; affective psychosis (bipolar disorder, or unipolar disorder with 
psychotic features); substance use disorders. Subjects with some re- 
creational substance use not associated with substance use disorder 
were deemed eligible for treatment. 

The enrolled patients received a three-years comprehensive, tai- 
lored and flexible intervention package. Prescribed tailored interven- 
tions included individual psycho-educational and motivational sessions, 
cognitive-behavioral psychotherapy, individual family psycho-educa- 
tion and support, therapeutic group activities (e.g., anxiety manage- 
ment, assertive and problem-solving training, etc.), various social group 
activities (e.g., music, multimedia, empowerment, computer training 
sessions, language classes, etc.), and supportive interventions on em- 
ployment, school, compliance with medication, and planning of re- 
creational activities (Cocchi et al., 2008; Meneghelli et al., 2010). 


2.2. Measures 


The following measures and indicators were used in this study: 
gender (boy or girl); age at first contact (years old); the 24-item Brief 
Psychiatric Rating Scale (BPRS; Overall and Gorham, 1962; 
Roncone et al., 1999); the DUP; past admissions to hospital for psy- 
chiatric reasons; recreational substance use not associated with sub- 
stance use disorder; family history of psychiatric disorders; dropout of 
treatment after the first two years for any reason; premorbid func- 
tioning. 

Gender and age at first contact were derived from a socio-demo- 
graphic form that is regularly filled in at first contact. The BPRS was 
used to assess general psychopathology, and it was the main measure of 
outcome. The BRPS is a 24-item measure of general psychopathology in 
a Likert format, with scores ranging from 1 (absent) to 7 (extremely 
severe). The range of possible BPRS total scores is 24 to 168, with 
higher scores indicating higher levels of psychopathology. The BPRS 
was regularly administered to the patients at inception and then every 6 
months, to assess levels of psychopathology and its change over time. 
BPRS was assessed by raters with a minimum of two years experience in 
rating patients diagnosed with psychosis. All raters were individually 
trained to apply the scales and the interviews that are part of the as- 
sessment package of the Programma2000. Inter-rater agreement among 
the clinical staff was regularly checked to assure good-to-acceptable 
concordance on the scales. Median intra-class correlation coefficients 
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(ICC) for BPRS, calculated with a 2-way mixed-effects model, were 0.70 
or above when performed on a small sample of patients (n = 25) during 
training. As in past studies (Ventura et al., 1993; Roncone et al., 1999), 
ICC were better for items related to patients’ self-report (> 0.80) than 
for items based on raters’ observation of the patients’ behavior (around 
0.60). Current guidelines suggest that ICC values between 0.5 and 0.75, 
between 0.75 and 0.9, and greater than 0.90 are indicative of moderate, 
good, and excellent reliability, respectively (Koo and Li, 2016). 

Improvement on the BPRS across time was measured as the differ- 
ence between the score at baseline (minus the minimum score of 24, as 
in Leucht, 2014) and the score at the end of the 2 years treatment 
(again, minus the minimum score of 24). 

The DUP was measured (in days) as the time elapsed from the onset 
of key symptoms (hallucinations, delusions or bizarre behavior) to the 
beginning of treatment (pharmacotherapy or psychotherapy) pre- 
scribed by a psychiatrist, on the basis of an interview that was ad- 
ministered to the patients at their fist contact with the Programma2000, 
the Early Recognition Inventory Retrospective Assessment of Symptoms 
checklist (ERIraos-CL), a 17-item screening checklist intended to select 
persons needing a more in-depth assessment. The ERIraos-CL has 17 
items designed to assist the exploration of individual proneness to 
schizophrenia, including items about changes in perception and 
thought interference and two items about paranoid ideation and hal- 
lucinations, which are more clearly indicative of psychosis 
(Hafner et al., 1992; Raballo et al., 2014). Data on past admissions to 
hospital for psychiatric reasons other than psychosis, on recreational 
substance use and substance use disorders, on family history of psy- 
chiatric disorders, and on premorbid functioning were based on a de- 
tailed interview with the patient and a key informant (a close relative, 
usually a parent). Past admissions to hospital for psychiatric reasons 
other than psychosis was assessed as yes/no and included any past 
admission to a psychiatric service as reported by the patients and 
confirmed by inspection of the records. Recreational substance use and 
substance use disorders were investigated by asking whether the patient 
had received any diagnosis of substance use disorder in the last 24 
months, or habitually consumed a drug from a list including cannabis, 
cocaine, heroin/opiates, amphetamine and its derivatives, hallucino- 
gens and a residual class of “others”. Recreational substance use was 
coded as yes/no whether the patient had a history of recreational 
substance use or not, irrespectively from the used drug. Family history 
of psychiatric disorders was investigated by asking if someone among 
close relatives (parents, grandparents, uncles and aunts, siblings, and 
cousins) had ever received a diagnosis of a mental disorder from a list 
including major mental disorders, or had ever been admitted to a 
psychiatric services, or had ever received the prescription of anti- 
psychotics, mood stabilizers or antidepressants as listed with their 
generic and commercial names. A family history of psychiatric dis- 
orders was coded as yes/no irrespectively from the number of affected 
relatives. Premorbid functioning was coded as yes/no according to 
whether the patients have had a decline in their functioning at school/ 
work or with their social relationship with peers in the past two years 
before contact with the center. Information on dropout was derived 
from the clinical charts and it was coded as yes/no irrespectively from 
the reasons for dropout. 


2.3. Statistics 


Exploratory data analysis and descriptive statistics were carried out 
using the Statistical Package for Social Sciences (SPSS) version 20. 
Additional analyses were carried out in R (R Core Team, 2018) using 
dedicated packages. All tests were two-tailed, with alpha set at p< 0.05. 
Means with standard deviations were reported for continuous variables. 
Counts and percentages were reported for categorical variables. Con- 
tinuous variables were tested with Student's t-test or ANOVA. Catego- 
rical analyses were carried out with the chi-square, along with Yates 
correction, whenever necessary. 
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A conditional growth model was used to analyze changes over time 
in the variable of interest (BPRS total score). The Restricted Maximum 
Likelihood estimation (REML), a linear mixed-effects model, was 
adopted using the “nlme” package running in R (Pinheiro et al., 2016). 
Change over time was tested in the variable of interest, as was the 
differential impact of gender (boys versus girls), age (less than 20 years 
old vs 2lyears old or older), DUP (less than 90 days vs more than 3 
months); past admissions to hospital for psychiatric reasons other than 
psychosis (yes vs no); recreational substance use (yes vs no); family 
history of psychiatric disorders (yes vs no); premorbid functioning 
(decline vs no decline); dropout of treatment after the first two years for 
any reason (yes vs no). 

The conditional (pseudo)R?, describing the proportion of variance 
explained by both fixed (time and group membership) and random 
factors (participants), was used as a rough measure of the model fit. The 
conditional (pseudo)R? was calculated according to 
Nakagawa et al. (2017). The analyses were performed with the 
“MuMIn” package running in R (Barton, 2018). 

Latent class growth analysis (LCGA) was applied to identify separate 
trajectories of subgroups of patients across time. LCGA identifies dis- 
tinct subgroups of individuals following a distinct pattern of change 
over time on the variable of interest (Andruff et al., 2009). According to 
the conventional growth model, intercept and slope of longitudinal data 
across time are expected to vary across individuals and random effects 
capture this heterogeneity. However, this approach assumes that all 
individuals are drawn from a single population with common para- 
meters. In LCGA, different groups of individual growth trajectories are 
expected to vary around different means (with the same or different 
forms), thus separate growth models are estimated for each latent class, 
each with its unique estimates of variances and covariate influences 
(Jung and Wickrama, 2008). In the model, the dependent variable 
(total score on the BPRS) was explained according to a quadratic 
function of time at the population level (fixed effects), and a linear 
function of time at the individual level. We expected that scores on the 
BPRS could decrease smoothly across time, or that scores could increase 
at the beginning and then decrease or, else, that scores could decrease 
at the beginning and then increase thereafter. One to six classes were 
tested. The package “lcmm” running in R was used to apply LCGA to the 
data (Proust-Lima et al., 2017). The best model was expected to mini- 
mize the values of the Akaike information criterion (AIC; 
Akaike, 1987), of the Bayesian information criterion (BIC; 
Schwarz, 1978), and of the sample-size adjusted BIC (SABIC, 
Sclove, 1987), and to maximize entropy, which was used to assess ac- 
curacy of participants’ classification (0 to 1), with higher values in- 
dicating better classification. Entropy values greater than 0.80 indicate 
a good separation of the identified groups (Ramaswamy et al., 1993). 
Participants were assigned to the latent class they had the highest 
probability of belonging to (average probabilities per class =>90% at a 
minimum was accepted). 

Logistic regression was used to assess the association between class 
membership and potential predictors of class membership. The same 
variables of interest that were tested in the conditional growth model 
were applied to the results of the LCGA. Differences between classes 
were expressed with odds ratio (95% confidence interval [CI]). 
Variance explained by the model (0 to 100%) was assessed with the 
pseudo-R? McFadden measure, with higher values being a reflection of 
better fit of the model (Long, 1997). 


3. Results 


Overall, there were 129 patients diagnosed with non-affective first- 
episode psychosis who completed the first two years of treatment at the 
Programma2000 during the interval of time that was considered in this 
study (1999 to 2016). Table 1 summarizes the main characteristics of 
the sample. Most patients were boys, and the mean age in the sample 
was 22 years old, with boys younger than girls (21 + 3 versus 23 + 4; 
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Table 1 
General characteristics of the sample (n = 129). 

All data are reported as mean (standard deviation); range, or counts (per- 
centage). 





Gender 
Boys 105 (81%) 
Girls 24 (19%) 


Age (years old) 22 (4); range: 16-31 


16 to 20 years old 53 (41%) 
21 years old or older 76 (59%) 
DUP (days) 152 (205); range: 0 - 730 
Less than 90 days 66 (51%) 
3 months or more 45 (35%) 
The DUP could not be determined 18 (14%) 
Past admissions for psychiatic reasons 
Yes 69 (53%) 
No 60 (47%) 
History of recreational substance use 
Yes 30 (23%) 
No 90 (70%) 
Not enough information 9 (7%) 
Family history of psychiatric disorders 
Yes 76 (59%) 
No 40 (41%) 
Not enough information 13 (10%) 
Decline in premorbid functioning 
Yes 95 (74%) 
No 23 (18%) 
Not enough information 11 (8%) 
Dropout of treatment after 2 years 
Yes 38 (29%) 
No 91 (71%) 
BPRS 
Baseline 50 (17); range: 24 - 105 
At 6 months 39 (14); range: 21 - 97 
At 12 months 38 (15); range: 24 - 103 
At 18 months 35 (9); range: 24 - 66 
At 24 monts 34 (9); range: 24 - 77 





t = 2.24, p = 0.026). About half of the sample had a DUP less than 90 
days, suggesting an acute onset of the condition. 

A past admission for psychiatric reasons other than psychosis was 
reported in 53% of the sample. Many patients reported a family history 
of psychiatric disorders (59%). Among those with a family history of 
psychiatric disorders, 24 reported that the diagnosis in the relative was 
a non-affective psychosis (32%), and 29 reported that the diagnosis in 
the relative was an affective disorder (38%), either major depressive 
disorder or bipolar disorder. Recreational substance use was observed 
in about a quarter of the sample (23%). Up to two-thirds of the parti- 
cipants showed evidence of a decline in premorbid functioning. Over 
time, the decrease in the BPRS score was on average 61%, suggesting 
robust improvement in psychopathology. 


3.1. Results of the conditional growth model 


Scores on the BPRS decreased over time, with a statistically sig- 
nificant decrease of overall scores on the BPRS across all tested models 
(Table 2). 

The decrease in BPRS scores was on average of 3 points every six 
months. Age, gender and DUP did not influence the time course of the 
decrease in BPRS scores. Those with past admissions for psychiatric 
reasons other than psychosis had a slower decline in BPRS scores 
compared to those without past admissions for psychiatric reasons. 
Patients who reported a family history of psychiatric disorders had a 
better course than those who did not. A history of recreational sub- 
stance use and the occurrence of decline in premorbid functioning be- 
fore the episode both were related to a lower improvement on the BPRS 
but the results were not statistically significant. Overall, the models had 
a good fit, with Conditional R? near or above 50%. Most of the variance 
was explained by the change over time of the scores on the BPRS and 
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the effect of the predictors was minimal. 
3.2. Results of the latent class growth analysis 


The indicators of fit decreased from the 1-class model to the 6-class 
model, however, the 4-class model had the best entropy, suggesting the 
best separation of the classes. The probability of assignment to their 
own class was above 90% in 86% of cases, which is the best perfor- 
mance among the models. Therefore the 4-class model was selected as 
the best description of the data in the sample (Table 3). 

In the 4-class model, one hundred and twelve participants were 
assigned to class 1 (87%), eight (6%) were assigned to class 2, five (4%) 
were assigned to class 3, and four (3%) were assigned to class 4 (Fig. 1). 

The majority of the patients showed a progressive decrease in the 
BPRS scores; a scant fraction showed a more stepped decrease from 
very high levels of psychopathology (class 2), or had a small im- 
provement at the first 6 months of treatment, followed by a worsening 
at 12 months and then a new decrease in BPRS scores (class 3). Finally, 
a residual group showed a sharp decrease in the first 6 months of 
treatment followed at 18 months by a worsening of their status (class 
4). 

No potential predictor was statistically related to the time course of 
BPRS scores as displayed across the four latent classes that were ex- 
tracted by the best LCGA model, essentially because of the small sample 
size of three classes (Table 4). 


4. Discussion 


Most patients that undergo treatment within an early intervention 
service are characterized by amelioration and ongoing changes. This 
result is congruent with the amelioration that was observed in patients 
with recent-onset psychosis after the start of antipsychotic treatment 
(Levine et al., 2010; Levine and Rabinowitz, 2010). Patients that have 
higher baseline scores of psychopathology have a more stepped de- 
crease of BPRS scores, probably requiring more intensive treatment. 
Overall, there was a trend for those with past admissions for psychiatric 
reasons other than psychosis, those with a history of recreational sub- 
stance use, and those who reported the occurrence of decline in pre- 
morbid functioning to show a slower decrease in BPRS scores. These 
patients might benefit from more intensive or more focused treatment 
to overcome their difficulties. Interestingly, those who reported a fa- 
mily history of psychiatric disorders had a better course than those who 
did not. This might depend on the family with experience of psycho- 
pathology being more sensitive to the impact of mental disorders, thus 
being more involved in the treatment of their kin when they have a 
psychotic breakdown. This is a topic deserving of further investigation. 

Overall, the improvement in the levels of psychopathology as 
measured by the BPRS was largely explained by the change over time of 
the BPRS scores, with minimal or null effect of the potential predictors. 
At least within the first 2 years of treatment, most predictors of the 
long-term course of schizophrenia-spectrum psychosis have no impact 
on the effects of treatment, or at least they have no impact on the pa- 
tients that are undergoing treatment within an early intervention ser- 
vice. Since the study didn't have a control group of patients receiving 
treatment as usual, we cannot firmly state that the treatment within an 
early intervention service removes most of the impact of the predictors 
of long-term outcome in schizophrenia. 

The Programma2000 is imbricated within the Lombardy mental 
health system of care. Lombardy is the largest and the most affluent 
region in Italy and has adopted an organizational model for the mental 
health care centered on multi-disciplinary teams, expected to be able to 
provide the full range of psychiatric care (Lora et al., 2012). In the last 
decade, a significant number of projects have been developed for early 
detection and interventions for severe mental disorders 
(Parabiaghi et al., 2019), and within these projects, Programma2000 
invested since the beginning into the improvement of patients’ 
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Table 2 
Results of the conditional growth models. 
Variables in the model Statistics 
Beta s.e. df t p Log likelihood (LL) Conditional R? 

Time =39 0.9 514 —4.4 0.0001 LL= — 2540.7, p = 0.020 44.3% 
Gender —0.7 4.3 127 —0.1 0.87 
Time x Gender 0.4 1.0 514 0.4 0.64 

Time —3.4 0.4 514 -7.1 0.0001 LL= — 2520.8, p = 0.012 49.6% 
Age 5.1 3.3 127 1.5 0.12 
Time x Age —0.9 0.7 514 -1.3 0.19 

Time —3.7 0.5 442 —7.4 0.0001 LL= — 2176.3, p = 0.007 47.9% 
DUP 5.5 3.5 109 1.6 0.11 
Time x DUP —0.9 0.8 442 -1.1 0.25 

Time -2.8 0.5 514 —5.0 0.0001 LL= — 2520.7, p = 0.0004 50.3% 
Past admissions 4.3 3.4 127 1.2 0.20 
Time x Past admissions RALO 0.7 514 -2.2 0.026 

Time =89 0.4 478 -9.5 0.0001 LL= — 2322.3, p = 0.024 52.0% 
Recreational substance use Time x Substance use 1.2 3.7 118 0.3 0.75 

0.1 0.8 478 0.1 0.08 

Time —2.6 0.6 462 -3.9 0.0001 LL= — 2260.9, p = 0.0001 49.8% 
Family history 6.1 3.6 114 1.7 0.095 
Time x Family history -1.9 0.8 462 —2.4 0.014 

Time -2.5 0.8 470 -3.1 0.0017 LL= — 2260.9, p = 0.0001 49.9% 
Prebordid functioning 9.7 4.1 116 2.3 0.019 
Time x preborbid functioning =1.6 0.9 470 -1.8 0.070 

Time =3.8 0.4 514 —8.4 0.0001 LL= — 2293.1, p = 0.0010 49.9% 
Dropout 2.1 3.7 127 0.5 0.57 
Time x Dropout 1.1 0.8 514 1.3 0.18 





Table 3 


Fit statistics for 1-6 class latent class growth mixture models. 











adherence to treatment, with a special emphasis, besides psycho- 
pharmacotherapy, on the regular attendance to psycho-educational and 
motivational sessions, cognitive-behavioral psychotherapy, individual 
family psycho-education and support, and therapeutic group activities 
(Cocchi et al., 2008; Meneghelli et al., 2010). Nevertheless, at Regional 
level, despite 72% of FEP patients receive a good multi-professional 
coverage, only one-sixth of them are involved in social and leisure ac- 
tivity, and activity devoted to family members are delivered to a min- 
ority of them and just 16% show adherence to the treatment regime 90 
days after the start (Lora et al., 2016). Poor adherence to psycho- 
pharmacotherapy during the first years after diagnosis influences an 
enhanced risk of hospital admission in FEP patients who are serviced 


4 Latent Classes 
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n. classes Log-likelihood AIC BIC SABIC entropy Posterior 
probabilities above 
90% in each class 
1 — 2514.88 5043 5063 5041 - - 
2 — 2495.13 5012 5043 5008 0.87 77% 
3 — 2478.54 4987 5029 4982 0.89 71% 
4 — 2471.56 4981 5035 4975 0.93 86% 
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Fig. 1. Treatment response trajectories over 2 years. On the left, raw data for all participants (each with their trajectory) and the estimated four-classes trajectories 
(wider lines). On the right, the smoothed trajectories of the four classes with the confidence interval, which is tighter in the larger class and larger in those with 
limited sample size. . 
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Table 4 
Baseline clinical variables across the 4 latent classes. 
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Class 1 Class 2 Class 3 Class 4 McFadden R? 
N = 112 N=8 N=5 N=4 
Gender 1.3% 
Boys 93 (83%) 6 (75%) 3 (60%) 3 (75%) 
Girls 19 (17%) 2 (25%) 2 (40%) 1 (25%) 
OR (95% CI) 1 1.6 (0.3 - 8.7) 3.2 (0.5 - 20.8) 1.6 (0.1 - 16.5) 
Age (years old) 0.3% 
16 to 20 years old 45 (40%) 4 (50%) 2 (40%) 2 (50%) 
21 years old or older 67 (60%) 4 (50%) 3 (60%) 2 (50%) 
OR (95% CI) 1 0.7 (0.1 - 2.8) 1.0 (0.1 - 6.3) 0.7 (0.1 - 4.9) 
DUP (days) 1.1% 
Less than 90 days 39 (41%) 4 (50%) 1 (20%) 1 (33%) 
3 months or more 56 (59%) 4 (50%) 4 (80%) 2 (67%) 
OR (95% CI) 1 0.7 (0.1 - 2.9) 2.8 (0.3 - 25.9) 1.4 (0.1 - 15.9) 
Past admissions for psychiatic reasons 5.9% 
Yes 58 (52%) 5 (62%) 5 (100%) 1 (25%) 
No 54 (48%) 3 (38%) 0 (0%) 3 (75%) 
OR (95% CI) 1 0.6 (0.1 - 2.8) 0 (0-0) 3.2 (0.3 - 31.9) 
History of recreational substance use 4.9% 
Yes 24 (23%) 4 (50%) 0 (0%) 2 (50%) 
No 79 (77%) 4 (50%) 5 (100%) 2 (50%) 
OR (95% CI) 1 0.3 (0.1 - 1.3) 0 (0-0) 0.3 (0.1 - 2.3) 
Family history of psychiatric disorders 1.5% 
Yes 
No 66 (66%) 5 (62%) 4 (80%) 1 (33%) 
OR (95% CI) 34 (34%) 3 (38%) 1 (20%) 2 (67%) 
1 1.1 (0.2 - 5.1) 0.5 (0.1 - 4.5) 3.9 (0.3 - 44.3) 
Decline in premorbid functioning 4.5% 
Yes 82 (81%) 7 (87%) 2 (40%) 4 (100%) 
No 19 (19%) 1 (13%) 3 (60%) 0 (0%) 
OR (95% CI) 1 0.6 (0.1 - 5.3) 6.4 (1.0 - 41.5) 0 (0-0) 
Dropout of treatment after 2 years 2.4% 
Yes 33 (30%) 1 (12%) 3 (60%) 1 (25%) 
No 79 (70%) 7 (88%) 2 (40%) 3 (75%) 
OR (95% CI) 1 2.9 (0.3 - 24.7) 0.3 (0.1 - 1.7) 1.2 (0.1 - 12.5) 





Latent Class 1 was used as a reference term. 


within the Lombardy mental health system of care; these patients rarely 
display proper use of the community service, in terms of regular at- 
tendance for receiving psychotherapy, psychoeducational activities, 
social supports, and multidisciplinary care (Corrao et al., 2015). Thus, 
even in the absence of a control group of patients receiving treatment- 
as-usual, it is reasonable to assume that the improvement observed in 
the sample is comparable to the superior improvement observed in past 
studies carried out in early intervention services, especially on positive 
symptoms (Correll et al., 2018). 

However, the improvement was not thoroughly and homogeneous. 
The LCGA has identified four classes in the treatment response trajec- 
tories of the sample, two of them with progressive improvement and 
two with some worsening over time. No potential predictor of outcome 
among those that were examined was specifically related to one of these 
classes, but the limited sample size might have prevented statistical 
estimation. Nevertheless, those patients incurring in a worsening of 
their status should be closely evaluated to identify factors that are re- 
lated to the worsening, including poor compliance to therapy, sub- 
stance use or the onset of a medical condition impacting on mental 
status. 


4.1. Strengths and limitations 


The use of state of the art statistics was the main strength of this 
study. Some limitations have to be acknowledged, too. The sample size 
was too limited to detect specific associations of the longitudinal classes 
with some of the predictors of outcome that were investigated in the 
study. Moreover, the assessment was executed at 6-month intervals, 
and this might have prevented some finer grounded definition of the 
trajectories. For some of the predictors, we had only binary data (yes or 
not), and this might have limited our chance of identifying relevant 


relationships. Another limitation is the lack of a control group. Indeed, 
it was not possible to agree with the psychiatric services operating in 
neighbor areas the application of an assessment detailed and repeated 
as the one that it is implemented in the Programma2000, thus no data 
on treatment-as-usual could be collected for comparison. 


5. Conclusions 


The study indicates that patients treated within an early interven- 
tion service are likely to undergo substantial improvement of their 
psychopathology. However, some of these patients may incur a wor- 
sening of their psychopathology during treatment and they need close 
reevaluation. 
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